A simple kinetic sensor to structural transitions 
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Driven non-equilibrium structural phase transformation has been probed using time varying re- 
sistance fluctuations or noise. We demonstrate that the non-Gaussian component (NGC) of noise 
obtained by evaluating the higher order statistics of fluctuations, serves as a simple kinetic detector 
of these phase transitions. Using the martensite transformation in free-standing wires of nickel- 
titanium binary alloys as a prototype, we observe clear deviations from the Gaussian background 
in the transformation zone, indicative of the long range correlations in the system as the phase 
transforms. The viability of non-Gaussian statistics as a robust probe to structural phase transition 
was also confirmed by comparing the results from differential scanning calorimetry measurements. 
We further studied the response of the NGC to the modifications in the microstructure on repeated 
thermal cycling, as well as the variations in the temperature drive rate, and explained the results 
using established simplistic models based on the different competing time scales. Our experiments 
(i) suggest an alternative method to estimate the transformation temperature scales with high ac- 
curacy, and (ii) establish a connection between the material-specific evolution of microstructure to 
the statistics of its linear response. Since the method depends on an in-built long-range correlation 
during transformation, it could be portable to other structural transitions, as well as to materials 
of different physical origin and size. 



I. INTRODUCTION 

Structural transitions constitute an evolving field of re- 
search comprising of unique and interesting mechanisms 
like creep and rupture^, plastic deformation^, externally 
induced crystal structure transformations*, the shape 
memory effect and pseudoelasticity 5 and so on. The rev- 
olution in material research, especially in nanoscale func- 
tional materials, ferromagnetic shape memory alloys etc., 
have added, both in complexity and diversity, to these 
transitions, particularly in the nature of their response 
to external stimuli such as mechanical stress, tempera- 
ture, current or magnetic field. 

A commonness in structural transitions in diverse ma- 
terials lies in the fact that the quenched disorder at the 
micro scale often manifests in discontinuities, avalanches 
or bursts and step changes in many of the physical re- 
sponses. A description of the response within the frame- 
work of continuum models, used for instance, to de- 
scribe the viscous flow at the onset of plasticity, breaks 
down, thereby demanding newer approaches to study 
these systems^. The crystallographic constraints and 
long range interactions lead to collective motions and 
avalanche-like behavior that have been the major focus of 
many recent works, indicating a universality in the statis- 
tical mechanics of these processes^—. They form a sub- 
set of a larger class of systems exhibiting non-equilibrium 
dynamic o 11 ^ 12 . These systems are found to posses a com- 
plex free energy landscape, where metastable states are 
separated by large potential energy barriers. 

Avalanches play a major role in driven non equilibrium 
systems manifesting in various forms like Barkhausen 
noise in magnetic hysteresis^, acoustic emission under 
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mechanical stress^ and so on. It has been observed that 
for such "crackling" noise, avalanches lead to fluctuations 
of physical properties that often become scale invariant 
or self similar. However, the statistics of these fluctu- 
ating dynamic quantities have largely been confined to 
the study of universality in the underlying physical pro- 
cesses, for example, in self organization of disorder—, but 
hardly exploited as a characterizing tool for material de- 
pendent transition properties themselves. Here, we aim 
to use higher order statistics of fluctuations as a sensor to 
structural transitions and the microstructural modifica- 
tions in a typical driven-nonequilibrium system, namely, 
the martensite transformation. 

Martensite transformation has been known to be a 
prototype of out-of-equilibrium phenomena 15 and at the 
same time constitutes the basis of a large variety of prac- 
tical applications including actuators, micro-pumps and 
so on in material systems, that display shape memory ef- 
fect and pseudoelasticity 5 . These transitions are known 
to be athermal, as large elastic energy barriers that sep- 
arate the metastable states require a driving force in the 
form of time varying stress or temperature for the tran- 
sition to proceed, and cannot be driven by thermal fluc- 
tuations alone—. Interplay of disorder and long range 
elastic forces is known to play a significant role in the 
transformation kinetics and it is this feature that we aim 
to exploit in the present study. 

The kinetics of martensitic transformation has been 
previously studied using acoustic emission measure- 
ments which probe the size distribution and duration of 
avalanches and their universal features&i&. However, the 
transformation temperature scales and the microstruc- 
tural evolution which are also essential for many de- 
vice applications are addressed only in limited number of 
cases^ 5 -. Generally, one resorts to the conventional meth- 
ods of calorimetry or standard resistivity measurements 
for the determination of these parameters 5 . Apart from 



2 




^50 200 250 ' 300 ' 350 



T(K) 

FIG. 1: (Color Online) (a) A typical resistance vs temper- 
ature plot exhibited by equiatomic Ni-Ti alloys. The tran- 
sition temperature scales are marked according to common 
practices, (b) Resistance vs temperature curves for the NiTi 
wire for various thermal cycles (tc). 

inherent inaccuracies and subjectivities involved, many 
of these techniques fail as one scales down the size and 
amount of materials^, becoming mostly inapplicable in 
the nanoscale systems. On the other hand, electrical 
noise, or fluctuations in resistivity, has also shown dis- 
tinct and interesting results in martensite transformation 
and avalanche dynamic s 18 ' 19 , and could be a potential 
tool to characterize the transition. Indeed, it was recently 
shown that higher order avalanche statistics is an excel- 
lent probe to martensite transformation, where detection 
of structural change depends on in-built long range cor- 
relation, rather than the amount of energy released^ 8 -. 
However, sensitivity of such a technique to small scale 
microstructural change is still not known. 

The present study is hence relevant in three ways. 
First, it provides a kinetic detector which could detect 
both magnitude and changes in the transformation tem- 
peratures as the material undergoes any slight modifica- 
tion in its structural morphology. Although resistivity 
measurements have been used for over three decades^ as 
a tool next to the calorimetric measurements, we show 
here that the phase transformation start and finish tem- 
peratures assigned based on the resistivity curves often 
lead to difficulties. We aim to establish that the fluctu- 
ation in resistance and its higher order statistics could 
reveal more detailed and accurate information about the 
transformation parameters as compared to conventional 
resistivity measurements. Second, though the athermal 
dynamics and the avalanche statistics have been well es- 
tablished in the context of crackling behavio r 8 ' 12 ' 20 , a 
clear link between this statistical approach and the mi- 
crostructural details of these materials is lacking. Statis- 



tics of fluctuations^ could again be used to probe this 
missing link, by altering the microstructure of the sam- 
ple through thermal cycling as well as by varying the 
drive rate. This is crucial to many applications of shape 
memory alloys since a minute change in the temperature 
scales on repeated use could result in major consequences 
in their applications. Third, even though the martensitic 
transformation is expected to be associated with a long- 
range elastic interaction, a direct probe to the nature and 
effect of this interaction is still absent. Here, we address 
this issue with higher-order fluctuation statistics as the 
experimental tool. A connection between the statistical 
parameters and the material specific microstructural de- 
tails that our results indicate might also be applicable to 
many forms of structural transition from a common plat- 
form. In addition, although the various timescales have 
been applied to explain the overlap of avalanches dur- 
ing the phase transformation, a clear link between these 
timescales and the material properties is not well estab- 
lished. We propose a simplistic model which accounts for 
the variation of the different transition parameters with 
increase in drive rate. 

Sections II and III of the paper summarize the material 
and the experimental methods respectively. In section 
IV, we present the major results and in section V, we 
present a brief discussion and conclude. 

II. MATERIAL 

The sample used was a nickel-titanium (NiTi) wire of 
rectangular cross-section 0.0005" x 0.002" procured from 
M/s Fortwayne Metals Inc. which was straight annealed 
and quenched in a strand furnace. The wire was found 
to show a two stage phase transformation from austen- 
ite (cubic B2) to R-phase (rhombohedral) and R-phase 
to martensite (B19/BI9' monoclinic) phase on cooling 
and vice versa on heating^. The wire was swept through 
a temperature range of 360 K-120 K for over 50 times, 
so that defect levels reached a stable configuration. Re- 
sistivity fluctuations or noise measurements were carried 
out at appropriate thermal cycles so that a correspon- 
dence could be drawn with the changes in microstruc- 
ture. The drive rate dependence was studied only after 
75 thermal cycles so that the effects of cycling and driv- 
ing induced changes in the sample morphology could be 
decoupled. 

III. EXPERIMENTAL METHODS 

A. Time averaged resistivity 

Resistivity measurements were carried out on the wire 
samples using a four probe geometry. A typical resistivity 
curve is as shown in Fig.l (a). The resistivity exhibits a 
hysteresis on cooling and heating indicative of the shape 
memory effect. The transformation temperatures can be 
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extracted using conventional formulation 5 as is shown. 
The wire exhibits a two stage phase transformation with 
an intermediate R-phase. The subscripts s and / denote 
the start and finish temperature scales of the different 
phases. 



B. Time dependent resistivity 

The sample was current biased using a low noise source 
and the voltage fluctuations were measured over an inter- 
val of time. We have used a five probe geometry and a 
dynamically balanced wheatstone's bridge arrangement 
(Fig. 2) to measure the noise as proposed by Scofield^i. 

Power spectral density (PSD) 

The power spectrum S(f) quantifies the noise in the 
frequency domain. If p(t) is the fluctuating resistivity, 
the power spectrum is defined as, 

w=T^(^)(£ dt ^w c " ar/t ) (i) 

Large classes of fluctuations in condensed matter sys- 
tems show a power spectrum that has 1/P dependence, 
with awl and hence called 1/f noise. 

In a typical experiment, the voltage fluctuations are 
acquired at a sampling rate of 1000 Hz. The data is sub- 
sequently decimated using a three stage decimation algo- 
rithm with appropriate anti-aliasing at each stage, which 
removes the unwanted signals. Care has been taken to re- 
move the extraneous noise in the system, using appropri- 
ate grounding and shielding using a faraday cage. An in 
depth review of the experimental setup and the measure- 
ment scheme including the basic digital signal processing 
steps could be found in our earlier paper—. 

Second spectrum and the non- Gaussian compo- 
nent 

Along with the power spectral density, we have also 
evaluated the higher order statistics of fluctuations in 
the form of second spectrum. The second spectrum is a 
fourth order moment of the fluctuations, that detects the 
weak correlations in the noise generating mechanism. We 
have used the approach followed by Seidler and Solin 22 , 
which is in turn a modification of the approach followed 
by Restle et ali2£. The resistivity fluctuations over a pe- 
riod of time, Ap(t), are recorded and then digitally band- 
pass filtered over a frequency band (/l, /jj). It is then 
squared point by point and the power spectrum of the 
time dependent variance of the noise signal is computed 
which is called the second spectrum, expressed analyti- 
cally as 

S (2) (/)=/ (Ap 2 (t)Ap 2 (t + T)) t cos(27rfr)dT (2) 
Jo 

The second spectrum was computed for each of these 
data sets, over a frequency band of fi = 1Hz and 
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FIG. 2: (Color Online) Schematic of the dynamically bal- 
anced wheatstone's bridge arrangement. Rl and R2 consti- 
tute the current limiting and balancing resistors respectively 
and rl and r2 constitute the sample with the middle con- 
tact grounded. The voltage is measured across a transformer 
preamplifier and then connected to a lock-in amplifier. 

fli = 3Hz, and was normalized with the square of the 
corresponding power spectrum, S p (f), as 

s {2 Hf) = S^(f)/l[ fH S p (f)df} 2 (3) 
J ft 

The Gaussian background has been calculated using the 
expression evaluated in Ref.22: 

(a) 2 1n(/ g -/)(/ L + /)// g / L 

Even weak correlations in the stochastic process con- 
tributes a non-Gaussian component (NGC) to the fluctu- 
ations, which can be readily identified from the deviation 
from its expected Gaussian background. In some cases, 
the NGC exhibits frequency dependence, as is clearly 
shown in Ref.18. Finite detection bandwidth sets the 
limit of sensitivity of s^(f) to non-Gaussian fluctua- 
tions. For comparison we calculated the total fourth or- 
der noise as, 

<? {2) = / ~ L S {2) (f)df (5) 
Jo 

Within a detection bandwidth of = 1 Hz and 
fli = 3 Hz, assuming S p (f) « 1//, the value of 
for a Gaussian background should be around 0.93 in ac- 
cordance with equation(4). In our case, we see that when 
correlations are not significant, we obtain a value close 
to 0.93 as the noise is approximately 1/f in nature over 
most of the temperature region. Any deviation from this 
Gaussian background can be considered to be the NGC, 
which is the signature of long-range correlations in the 
system. 
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IV. RESULTS 
A. Resistivity fluctuations and statistics 

We first probed the signature of avalanches in the resis- 
tivity fluctuations across the transformation zone (using 
temperature or stress in the case of martensite transfor- 
mation). Since the martensite transformation in NiTi 
is primarily athermal, the avalanches are induced only 
by a continuous sweeping of T. The characteristics of 
transformation, and hence the statistics of fluctuations, 
depend crucially on the sweep rate r = dT/dt, and will 
be discussed in detail later. Fig. 3(a) shows the time- 
dependence of resistivity around T = 250 K for a NiTi- 
device cooled at r = 1 K/min. It is clear that the fluc- 
tuations are much smaller in the static condition, i.e. 
keeping T constant at the same value, when avalanches 
are expected to be absent. 

In a typical experiment, the fluctuations in p were mea- 
sured with T being swept from T Af ~ 325 K to 
T -C Mf ~ 180 K (Af and Mf are the Austenite finish 
and Martensite finish temperature scales respectively), 
with the T-range divided in equal time windows of 8.74 
minutes. Both PSD and higher order statistics of noise 
were subsequently calculated in each window and plotted 
against the temperature corresponding to its midpoint. 
The PSD of the fluctuations varied as S p (f) ~ 1/ f a over 
the experimental bandwidth, where a ~ 1.1 — 1.3 (See 
Fig. 3b). The magnitude of a was found to be slightly 
smaller than that in NiTi thin films observed earlier, 
which could be due to greater inhomogeneity in the bulk 
NiTi system in the present case. The PSD magnitude 
was found to be T-dependent as well, increasing mono- 
tonically with T in all devices over the observed T range 
(See Figs.6(a) and 8(a)). 

The behavior of the second spectrum was found to be 
different. Typically, S*- 2 * 1 agrees with the Gaussian back- 
ground (the dashed line in Fig. 3c), except at certain 
values of T where it was found to be orders of magnitude 
higher, indicating significant non-Gaussian contribution 
to noise. Fig. 3c illustrates this for a few representative 
T, where s 1 - 2 ' at various temperatures T < 240A" is not 
only much higher than the Gaussian background, but 
also develops a weak slope ~ 1/ f 13 (J3 ~ 0.4 — 0.6). 
The absolute magnitude of was however found to be 
highly sample dependent, without any clear trend, but its 
variation with T was qualitatively similar in all devices, 
which we shall now focus on. No non-Gaussian compo- 
nent of noise could be detected in the static mode (not 
shown), implying non-Gaussian behavior to arise from 
correlated avalanches (not shown). 

Unlike PSD, we find the NGC to be non-monotonic in 
T, and manifest as strong peaks in at specific temper- 
atures. One expects four structural transitions in NiTi 
within this range: Austenite to R-phase (R«— A) and R- 
phase to martensite (M<— R) during cooling, and marten- 
site to R*-phase (M->R*) and R*-phase to austenite 
(R*— >A) during heating. In Fig. 4(a) (cooling) and 4(c) 



(heating) these four transitions are identified by the la- 
tent heat release in the differential scanning calorimetry 
measurements (Fig. 4(b) and 4(d)). Strikingly, the peaks 
in cr( 2 ) appear primarily around these phase transforma- 
tion zones. Similar correlation of non-Gaussian electri- 
cal noise and first order structural phase change was ob- 
served in NiTi thin film as welli^, outlining an excellent 
kinetic detector of structural phase transition. The NGC 
in resistivity noise has been attributed to emergence of a 
long range correlation in many systems with electronic 2 ^, 
structural^ or magnetic 25 phase transitions. The magni- 
tude of cr^ 2 ) is directly associated to a correlation length 
£ which, in case of martensitic alloys, characterizes the 
range of elastic correlation in the system. As the transi- 
tion is approached, the avalanches become more frequent, 
which overlap in space, as well as in time in the case of 
faster sweeps, leading to strongly hierarchical and coop- 
erative atomic kinetics. The non-Gaussian fluctuations in 
p is a result of electron wave function coupling to the size 
of atomic avalanche, i.e. the number of atoms undergoing 
displacive movement. Consequently, £ may attain a max- 
imum at a critical temperature (T c ) where the avalanche 
overlap is also most probable. This would occur when 
maximum amount of parent phase is transformed into 
the product phase, as confirmed by the agreement of T c 
with the peak heat release/absorb temperature scales in 
DSC measurements. 

In the following, we stress on two specific questions: 
(1) Can NGC provide a convenient scheme to evaluate 
the transformation temperatures, such as the Mf, Af 
etc? Secondly (2) how do small scale microstructural 
modifications influence the non-Gaussian characteristics? 
The microstructural evolution of the material depends 
mainly on two processes: Thermal cycling, which changes 
the local disorder level, and also the drive rate, which 
leads to a significant modification in the internal stress 
levels and avalanche overlap. To analyze these cases in 
detail, we shall focus mostly on a specific transformation 
zone, namely that from R-phase to the martensite phase 
(M-<— R) while cooling, which in most devices extend from 
T ~ 180 K<- 240 K. 



B. Transition temperatures 

From the perspective of material characterization, re- 
sistance measurements have long been used as a conve- 
nient probe to measure the phase transformation tem- 
perature scales. In case of martensites, these include M s 
(martensite start), Mf (martensite finish), A s (austen- 
ite start), Af (austenite finish) etc. However, as can be 
seen from Fig. 5, where p — T data from three NiTi wires 
are compared, it often becomes very difficult to assign an 
unambiguous start and finish temperature to the phase 
transformation region. Instead, by comparing the range 
of T in DSC data over which the avalanches are observed, 
to the width of the peak, an alternative scheme can 
be suggested. As can be seen in the context of marten- 



5 




0.02 



FIG. 3: (Color Online) (a) A typical time series of the resis- 
tance fluctuations while the NiTi sample was being ramped at 
lK/min around a mean temperature of 250 K and while keep- 
ing the temperature constant at 250 K (b) A typical power 
spectrum showing the 1//" behavior, (c) The normalized 
second spectra for various temperatures during a cooling run. 
The increase in the integrated second spectra, which ulti- 
mately leads to a peak in <j' 2 ' is seen as a shift in the spectra, 
in this plot. The Gaussian background is plotted as a dotted 
line. 




FIG. 5: (Color Online) Integrated value of the normalized 
second spectra a*- 2 ' as a function of temperature for three dif- 
ferent samples shown along with the corresponding resistance 
plots. The heating and cooling directions are shown using 
appropriate arrows. The region over which a^ 2 ' deviates from 
the Gaussian background (dotted line) is marked using ver- 
tical lines. The point at which cr' 2 ' reaches its maximum is 
defined as the critical temperature T c . 



site transition (M-s— R) in Fig. 5, the higher and lower T 
(denoted by the vertical dashed lines) at which the value 
of cr^ 2 ) starts to deviate from the Gaussian expectation 
could be assigned as M s and Mf respectively. In Fig. 5a 
the A s and Af are identified by the vertical dashed lines, 
locating the R<s— A transition. The same scheme can be 
adopted for the other transitions as well. 




T(K) 

FIG. 4: (Color Online) The integrated value of the normalized 
second spectra a* 2 ' as a function of temperature showing the 
peaks corresponding to the transitions for (a) cooling and (c) 
heating. The corresponding DSC plots are shown in the lower 
panels (b) and (d). The calculated Gaussian value integrated 
over the frequency bandpass is shown as a dotted line. The 
details about the transforming phases are denoted in (a) and 
(c). 



C. Thermal Cycling 

The system was swept through a temperature range of 
360 K-120 K. The corresponding resistance vs. temper- 
ature plots are as shown in Fig. 1(b). It is clear that on 
thermal cycling the defect levels increase and stabilize af- 
ter around 50 thermal cycles. Subsequently, we measured 
resistance fluctuations while the system was being driven 
at 0.15 K/min across the temperature range of 240 K- 
180 K (cooling). This region is where the transformation 
from the R-phase to martensite phase takes place. It was 
observed that the spectra exhibits l/f a behavior. The 
frequency exponent a for the entire temperature zone for 
the various cycles is plotted in Fig. 6(a), which is found to 
cluster around 1.2 and does not show any prominent vari- 
ation with cycling. Fig. 6 (a) also shows a clear decrease in 
the noise magnitude with decrease in temperature. The 
noise magnitude however shows no systematic cycle de- 
pendence, although some earlier studies indicated a weak 
increase in its magnitude with cycling. 

Strikingly however, we observe the critical point T c {= 
T^^ R ) in shows a clear shift towards lower tem- 
peratures on thermally cycling the sample. To make this 
behavior clearer, we have plotted T c as a function of ther- 
mal cycles in Fig. 7. To compare with the changes in the 
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FIG. 6: (Color Online) (a) Power spectral density (open sym- 
bols) and the frequency exponent a (closed symbols) as a 
function of temperature for various thermal cycles, (b) The 
variation of as a function of temperature for different 
thermal cycles. The Gaussian expectation is plotted as dot- 
ted lines. The deviation from the Gaussian behavior of the 
second spectrum around the critical temperature is evident. 
The fits shown are Gaussian. 



resistance with respect to cycling, we have also plotted 
the resistance value at 245 K during cooling for as a func- 
tion of cycles in the same plot. It is clear that both T c 
and resistivity stabilize after similar number of cycles, 
indicating a common underlying physics. 



D. Drive rate dependence 

To explore the drive rates effects, the system was 
driven at five dive rates r = 0.15 K/min, 0.3 K/min, 
0.5 K/min, 0.7K /min and 1 K/min. To decouple the 
effects of drive rates from thermal cycling, this set of ex- 
periments was done only after the disorder levels were 
stabilized to the equilibrium limit. It was observed that 
the peak position T c for 0.15K/min after and before the 
drive rate dependence tests were exactly the same, indi- 
cating that the drive rate experiments have not affected 
the system's morphology in an irreversible manner. The 
PSD docs not show rate dependence (Fig. 8(a)), which 
is essentially a feature of the athermal nature of the 
phase transition. The frequency exponent a is plotted 
in Fig. 8(a), which shows no rate dependence either. 

Fig. 8(c) shows the variation of for the various 
rates. As the drive rate increases, there is a gradual shift 
of T c towards lower temperatures. As the drive rate be- 
comes larger, the number of sampling windows decreases, 
which lead to a larger error bar. In addition, the magni- 
tude of is found to increase on increasing the rate. 
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FIG. 7: (Color Online) The variation of T c as a function of 
thermal cycles calculated from the Gaussian fits shown in 
Fig. 6(b). The curve is a guide to the eye. The resistance 
value at 245K during cooling for the various thermal cycles is 
also plotted for comparison. The inset shows the variation of 
cr' 2 ' at T c for the various cycles. 
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FIG. 8: (Color Online) (a) Power spectral density (open sym- 
bols) and the frequency exponent a (closed symbols) as a 
function of temperature for various drive rates, (b) The vari- 
ation of <t' 2 ) as a function of temperature for the various 
drive rates. The Gaussian expectation is plotted as dotted 
lines. The deviation from the Gaussian behavior of the sec- 
ond spectrum around the critical temperature is evident. The 
fits shown are Gaussian. 



These two features are shown in Fig. 9. 



V. DISCUSSION 

The experimental data indicates that both thermal cy- 
cling and drive rates have very similar effects on the be- 
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FIG. 9: (Color Online) The critical point T c calculated from 
the Gaussian fits in Fig. 8(b), as a function of drive rate. The 
magnitude of cr' 2 ' at T c as a function of drive rate is plotted 
on the right hand side axis. The curves are guide to the eye. 



havior of the resistance fluctuations and its higher order 
spectra. In this section, we would work to develop a 
qualitative understanding of these effects based on the 
evolution of elastic energy and overlap of avalanches as 
the system traverses through thermal cycles and the vari- 
ations in drive rates. 
Thermal Cycling 

The effects of thermal cycling have been extensively 
investigated by Miyazaki et al. in NiTi samples^. It has 
been established that on cycling dislocations are intro- 
duced in samples which are not extensively heat treated 
prior to cycling. This generic behavior has been charac- 
terized using resistance measurements (Fig. 1(b)), which 
shows an increase in resistivity of the wires on thermal 
cycling due to the increased number of dislocations. How- 
ever, the initial increase in dislocations causes an inter- 
nal strain field which then opposes further proliferation 
of dislocations leading to a saturation of disorder level. 
Perez-Reche et al. have attributed this feature to a self- 
organization of defects^ 4 -. 

To explain the shift in T c , it is essential to understand 
the changes in the internal elastic strain energy and its 
implication on the transition temperature scales. From 
the DSC measurements as well, it has been observed that 
the martensite start temperature M s decreases on ther- 
mal cycling. This has been attributed to the increase 
in dislocation stress fiel d 26 ' 27 , which results in increas- 
ing the potential energy barrier between the metastable 
states and hence require an additional under-cooling to 
drive the phase transformation. This scenario readily ex- 
plains the variation in T c as well, but the absolute magni- 
tude of cr*- 2 ) was found to be insensitive to thermal cycling 
(inset of Fig. 7). This implies that while the local elas- 
tic energy profiles are modified during thermal cycling, 
the extent of long range correlation length remain largely 



unaffected. 

Drive Rate Effects 

The effect of drive rates has been probed in our earlier 
work as well 18 and Fig. 8(a) shows very similar results. To 
understand the second spectral characteristics (i.e. the 
magnitude of and the deviation in T c ), we adopt a 
formalism similar to the one followed in Ref.28. 

A driven structural phase transition can be described 
using three competing timescales: (1) The avalanche re- 
laxation r av , which depends on the interaction in the 
system and the disorder level, being typically ~ 10~ 8 — 
10 -4 sec in most systems. (2) The time scale of thermal 
fluctuations (Tth), which depends on the elastic energy 
barriers separating the metastable states, and depends 
on the mechanisms for the transition like slipping, twin- 
ning etc. In an athermal transition the effects of thermal 
fluctuations are small, and r t h is generally taken to be 
very large (r t h > 10 5 — 10 6 sec) in conventional shape 
memory alloys such at nitinol. (3) The drive rate r itself 
defines a time scale Tdr <x 1/r, which is well separated 
from both Tth and r av in most cases, which quantifies the 
athermal limit as Tth/Tdr — > oo, and the adiabatic limit 
as T av /Td r — > 0, where the finite drive does not affect the 
characteristics and propagation of avalanches. 

In a real system both Tth and T av are finite, and at finite 
drive two processes can occur: (i) the activation of the 
transition by thermal fluctuations, and (ii) the overlap of 
avalanches. The behavior of both and T c shown in 
Fig. 9 can be explained with these two processes, assum- 
ing in general ry, > Td r > T av and a scenario depicted 
schematically in Fig. 10. The vertical dashed line repre- 
sents a particular sample with material specific Tth and 
r av . As drive rate increases, r av /Td r also increases, re- 
sulting in larger probability of avalanche overlap. Due 
to the long range nature of elastic force, superposition 
of avalanches would couple events separated in time and 
space, leading to longer correlations, and hence increase 
in the magnitude of as seen in Fig. 9. 

The variation in T c , on the other hand, can be as- 
sociated to an increasing Tth/Tdr when r increases, i.e. 
as the system is effectively driven towards the athermal 
limit. Hence the number of smaller avalanches triggered 
by thermal fluctuations decreases, requiring an additional 
undercooling to transform the same amount of the par- 
ent phase to product. This naturally leads to a lower 
value of T c . It is to be noted that the effect of finite drive 
rate on T c would be insignificant if the system is truly 
athermal Tth/Tdr oo, which is traditionally assumed 
for the martensite transformation in NiTi 2 ^. However, 
our experiments here indicate that the thermal fluctua- 
tions do play a minor role in NiTi. It is also important 
to note that the power spectral density of noise changes 
very little for various drive rate a 18 ' 19 , indicating the sec- 
ond spectral properties to be far more sensitive to small 
deviations from athermal behavior and adiabatic limit 
than the noise itself. 
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^ Athermal 




FIG. 10: (Color Online) Schematic representation of the sim- 
plistic model based on the competing timescales, described 
in the discussion section. r\ , r-z and rz denote three different 
rates. 



VI. CONCLUSION 

In this work we have adopted a novel method based 
on higher order statistics of fluctuations to probe into 
the kinetics of martensite structural transformation in 
nickel-titanium alloys. The non-Gaussian component of 
noise in electrical resistivity of these systems is found to 
be an excellent probe, which relies on the built-in long 
range elastic correlations in this system. The sensitivity 
of the technique to small scale microstructural modifica- 
tions was also addressed by measuring noise as a func- 
tion of thermal cycling and the drive rate. We believe this 
could serve as an independent and sensitive tool to probe 
various forms of structural transition in a wide variety of 
(electrically conducting) material systems. 
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